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ABSTRACT 

We show, via a simple parametrization of the multiplicity distribution of 
charged particles in e + e~ annihilation at the Z° peak in terms of the weighted 
superposition of two negative binomial distributions, that both the shoulder 
structure in the intermediate multiplicity range and the oscillation in sign of 
the ratio of factorial cumulants over factorial moments of increasing order are 
related to hard gluon radiation. 
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We show, via a simple parametrization of the multiplicity distribution of charged 
particles in e + e — annihilation at the Z° peak in terms of the weighted superposi- 
tion of two negative binomial distributions, that both the shoulder structure in the 
intermediate multiplicity range and the oscillation in sign of the ratio of factorial 
cumulants over factorial moments of increasing order are related to hard gluon 
radiation. 

1 Common origin of shoulder structure and moments' oscillations 

One of the main still open problems in multiparticle dynamics is that of attain- 
ing an integrated description of final particle Multiplicity Distributions (MD's) 
and of the corresponding correlation functions properties. One expects that 
features detected in terms of one of the two observables have a sound physical 
counterpart in terms of the other. It is of course a quite difficult task to ex- 
plain facts occurring in the two domains by means of the same physical cause 
and accordingly to show that they have a common origin. We discuss in the 
following a successful example of this search in e + e~ annihilation. Here two 
interesting features are observed at the_Z2. energy: the MD shows a shoulder 
in the intermediate multiplicity rangelira'tm and the ratio of factorial moments 
to factorial cumulant moments changes sign as a function of its orderoO The 
two observables are strictly linked, as factorial moments, F q , and factorial cu- 
mulant moments, K q , can be obtained in general from the MD, P n , through 
the relations: 
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Figure 1: The ratio of factorial cumulant moments over factorial rwciments, H q as a function 
of q; experimental data (diamonds) from the SLD Collaboration □ are compared with the 
predictions of several parameterizations, with parameters fitted to the data on MD's: a full 
NBD (dotted line); a truncated NBD (dot-dashed line); sum of twpfull NBD's as per Eq. ^ 
(dashed line); sum of two truncated NBD's as per Eq. N (solid line). 



and 

^ = ^-E( i A V^- (2) 

The ratio of the two mentioned quantities: 

H q = K q /F q (3) 

is well suited to theoretical and experimental studies!] 

QCD calculations are available for both the MDLj and the ratio H q u Both 
fail to reproduce the data after a simple application of Local Paoton Hadron 
DualityH as hadronization prescription (or its generalized versionp to be used 
with moments of higher orderE3). The most recent calculation of the MD has 
been done in MLLA gluodynamicsQ KNO scaling is present only asymptoti- 
cally and the predicted shape at current energies is much narrower than in 
previous calculations!^ but still the head and the tail of the distribution can- 
not be correctly reproduced. The ratio H q has also been calculated beyond 
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DLAp and has been predicted to show a negative minimum around q = 5 and 
then to oscillate in sign as a function of the order q; this agrees qualitatively 
with the behaviour seen in the datatj 

Many phenomenological parameterizations have been used to describe data 
on MD's, the most common (and simplest) being the Negative Binomial Dis- 
tribution (NBD)E3 and the Log-Normal Distribution (LND)N Both work well 
at lower energies, but fail to reproduce the shoulder structure at y/s = 91 GeV. 
In particular this failure can be seen if the residuals (difference between the 
data and the parametrization divided by the error on the data) are examined: 
a clear structure is seenP The SLD Collaboration hasi-also shown that these 
parameterizations fail in the sector of. the H q momentsp even after taking into 
account the effect of finite statistics til This is particularly evident in Figure [y, 
if one looks at the dotted (full NBD) and dot-dashed (truncated NBD) lines. 

The shoulder structure was first explained by the Delphi Collaboration: lil it 
results from the superposition of samples of events with a fixed number of jets. 
In each sample taken separately there is no shoulder but, since they have a 
different average multiplicity, a shoulder appears in the full sample. In addition 
it was found that the NBD .describes well the MD in these samples of events 
with a fixed number of jetsEj It should also be remembered that a shoulder 
has been seen in pp collisions at yfs = 90Q-CeV by the UA5 Collaboration^ 
and it has been confirmed at the TevatrontZl A good fit to the UA5 data was 
performed with the weighted sum of two NBD'slia 

Following the above observations, we propose a parametrization of the MD 
which is the weighted sum of two components, one to be associated with 2-jet 
events and one to be associated with events with 3 or more jets. The weight in 
this superposition is then the fraction of 2-jet events, which is experimentally 
determined, it is not a fitted parameter. This decomposition depends of course 
on the definition of jet; in particular it depends on the particular jet-finding 
algorithm, and on the value of the parameter y m i n which controls the algorithm. 
The Delphi Collaboration has used the JADE algorithm and published values 
for the 2-jet fraction and for the MD's at y min = 0.02, 0.04, 0.06, 0.08. We 
have performed the fit for each of these values. 

Concerning the experimental data, it should be noticed that the experi- 
mental errors on each point of the published MD's are correlated with adjacent 
bins: in our fit we cannot take this correlation into account, therefore it is not 
possible to compare directly the values of the \ 2 we obtain with those obtained 
by the experimental collaborations. Furthermore, the extraction of H q from 
the published P n also suffers from a similar problem, and the. errors we show 
in the figures are estimates obtained by a statistical method t3 (except for the 
data of the SLD Collaborationta). 
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As for the particular form of the MD in the two components of our fit, we 
have chosen the NBD, because it has successfully been fitted to the data for 
the samples of events with fixed number of jets. In practice we perform a fit 
to the MD's with a four parameter formula: 



Here P„ BD (n, k) is the standard NBD of parameters n and fc; notice that we 
have taken into account the charge conservation law, which requires the final 
charged particle multiplicity to be even. The proportionality factor is fixed by 
requiring the proper normalization for P n . 

Results of our fit to the data of four experiments are shown in Table [l| 
we find \ 2 P er degree of freedom equal to or smaller than 1, and values of the 
parameters consistent between different experiments; they are also consistent 
with those obtained by the Delphi Collaboration in fitting their 2-jet and 3- 
jet data separately with a NBD. These findings are visually summarized for 
a = 0.767 in Figure ^|, where the residuals do not show structures. 

In Figure |3j we compare the experimental data on HqS with the values 
obtained from a formula that takes the truncation effect into account, too: 



where n m i n and n max are the minimum and maximum observed multiplicity, 
and a proportionality factor ensures proper normalization. We obtain a very 
good agreement with the data. Notice that it is not possible to reproduce the 
behaviour of the ratio H q without taking into account the limits of the range 
of the available data. This can be seen in Figure [l], by comparing the dashed 
(full formula) and the solid (truncated formula) lines. 

In conclusion, the observed behavior of H^s results from the convolution of 
two different effects, a statistical one, i.e., the truncation of the tail due to the 
finite statistics of data samples, and a physical one, i.e., the superposition of 
two components. The two components can be related to 2- and 3-jet events, i.e., 
to the emission of hard gluon radiation in the early stages of the perturbative 
evolution. 
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Table 1: Parameters and \ 2 pe»-iDegree ofj-iFreedom (df) of the fit to experimental data 
from ALEPHg DELPHltJ OPALEI and SLDq Collaborations with the weighted sum of two 
NBD's. Results are shown for different values of a corresponding to the fraction of 2-jet 
events, /, experimentally measured by DELPHI Collaboration |^£| at different values of the 
jet-finder parameter y mln . NBD parameters extracted by the DELPHI Collaboration by 
fitting MD's of samples of events with 2- and 3-jets at different values of y mln are also shown 
for comparison in the last column. 
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Figure 2:n Chaffed particles-iMD's in full phase space, P„, at the Zn peak from ALEPH|j 
DELPHltl SLD □ and OPALB Collaborations are compared with Eq. | with a = 0.767 (see 
Table h for the values of the corresponding parameters) (solid lines); dotted lines indicate 
the two separate NBD contributions. The lower part of each plot shows the residuals, R n , 
i.e., the difference between data and theoretical predictions, in units of standard deviations. 
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Figure 3: The ratio of factorial cumulant moments over factorial moments, H q as a function of 
q; experimental data (diamonds) from ALEPH, DELPHI, SLD and OPAL Collaborations are 
compared with Eq. fej for different values of a (see Table 1 for the values-pf the corresponding 
parameters). In the figure only statistical errors of SLD dataQ are shown. 
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